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Abstract 

We consider two aspects of the ground state of the U = oo Hubbard ladder: 
ferromagnetism and the metal-insulator transition at quarter-filling. First, we 
present rigorous results for the U = oo Hubbard ladder in the limit of the large 
inter-chain hopping (t±/t — * oo). In this limit, the total spin S of the ground 
state is shown to be zero for the electron density n < 0.5 and its maximum 
(S = 5 max ) for n > 0.5. The charge gap at quarter-filling is 2t±. We extend 
these results to finite t±/t by means of the density- matrix renormalization 
group method. We estimate the phase boundaries with respect to spontaneous 
magnetization and the charge gap at quarter-filling for finite t±/t. Applying 
the extended Aharonov-Bohm method, we give numerical evidence that the 
critical ratio t±/t, above which the charge gap opens, is less than 0.001. 
Ferromagnetism in the t-J ladder is briefly discussed. 

PACS numbers: 71.27.+a, 75.10.Lp, 71.30.+h, 73.61.At 



1 



Typeset using REVTfjX 



I. INTRODUCTION 



Recently, a lot of attention has been given to the Hubbard or t-J ladders. One of the 
reasons is that these models are believed to describe the essential features of various materials 
such as (VO) 2 P 2 7 [0 or SrCu 2 3 [@. Another reason is that some properties of the ladder 
models are closer to planar models than to single chain models. One such property is the 
enhancement of the d x 2_ y 2-\ike pairing correlation. In relation to high-T c superconductivity 
HQ , this feature and the existence of the spin gap have been studied by various approaches 
such as bosonization [0-0], projector Monte Carlo |§, exact diagonalization fl-pl and the 
density-matrix renormalization group (DMRG) method [T^jT^ . 



Another interesting aspect of the ladder models is the existence of the ferromagnetic 
ground state. Nagaoka's theorem |T3|-|nj holds for the ladders and the planar models, 
although not for the single chains. Hence, ferromagnetic ground states may exist in the 
strong-coupling regime near half-filling for the ladders and the planar models, but not for 
the single chains. For two dimensions, many workers have investigated the ferromagnetic 
ground state at U = oo for finite hole-density |T7H24| , in order to clarify the origin of itinerant 
ferromagnetism from the electron-correlation view-point. In spite of such attempts, results 
are inconclusive, because of the lack of efficient methods. On the other hand, for the ladders, 
there is an efficient method, i.e. the density-matrix renormalization group (DMRG) method 
proposed by White p5|] . Liang and Pang |24]] applied this method to the U = oo Hubbard 



ladder for t±/t = 1, and obtained some indications of ferromagnetism. As for the (two- leg) 



ladder, the DMRG calculation in ref. [24] suggests that the fully-polarized ferromagnetic 
state is one of the ground states for 5 < 0.22 (5 : hole density) and that the ground state 
is a spin-singlet for 5 ~ 0.4. They have tried to extend these results to two dimensions by 
investigating multi-leg ladders. In this paper, we extend their results to various values of t± 
and J for the t-J ladder, in order to understand the ground-state properties of the Hubbard 
or t-J ladders in the strong-coupling regime. 

Another interesting feature of strongly correlated electron systems is the metal-insulator 
(MI) transition. For single chains and planar models, there have been a lot of works on 



the MI transition fl26|-pT|j. However, for ladder models, there are relatively few. One of the 
characteristic features of ladder models is the band structure. In the weak-coupling regime 
(U ~ 0), the low-energy physics for n < 1 is effectively described by the bonding band, if 
t±/t is large enough. Thus, the MI transition in this parameter regime is essentially the 
same as that of the single chains 0. This argument may be true for t± ^> U. A natural 
question arising here is what happens in the opposite limit, i.e. U 3> t±. In order to answer 
this question, we consider the MI transition in the U = oo Hubbard ladder for simplicity. 

In the present paper, we mainly consider the U = oo Hubbard ladder (or t-ladder), 
which is defined as the t-J ladder at J = J± = 0. Later, in Sec.[TTI|, we briefly discuss 
ferromagnetism of the t-J ladder. The Hamiltonian of the t-J ladder is defined as follows : 

H t j = H t + Hj, (1.1) 

n t = -tJ2(4M +la + ^c 2 i+la + h.c.) 

ia 

-t±£(c£c£ + ii.c.), 
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Hj = JJ2(S} ■ Sl +1 - -n}n] +l + S> • S* +1 - 

i 

+J±Il(Sl'S>-±nln*), 

i 

where c"} denotes a creation operator of an electron at site i with spin a (a =j, j) in the a-th 
chain (a = 1, 2) with the constraint that no site is doubly occupied, i.e. c^J = c"J(l — n"_ a ). 
The number operator n^. is defined as nf a = c*}cf a , using the standard electron creation 
operator c*}. The spin operator at site i in the a-th chain is defined as = \ YL^ ( ^l (7 ^ c %i 
where er is the vector of Pauli matrices. The number of electrons, the number of sites and 
the number of rungs are denoted by N e , N s and L(= N s /2), respectively. The electron 
density n and the hole density 5 are defined as n = N e /N s and 8 = 1 — n, respectively. The 
maximum value of the total spin S is denoted by S' max (= N e /2). 

This paper is organized as follows: In Sec [IT], we present rigorous results on the ground 



state of the t-ladder in the limit t±/t — > oo. In Sec PH we present numerical results on 



ferro magnetism in the t-ladder and the t-J ladder. In Sec [TV], we discuss the metal- insulator 



transition at quarter-filling for the t-ladder. Section is a summary. 



II. T-L ADDER IN THE LIMIT T ± /T -> oo 
In this section, we prove the following statements: 

1. The ground state of the t-ladder in the limit t±/t — > oo for n < 0.5 is a spin-singlet 
(S = 0) and is unique in finite-size clusters with an even number of electrons with 
open boundary conditions. 

2. The ground state of the t-ladder in the limit t±/t — > oo for n > 0.5 has the maximum 
total spin (S = S max ) and is unique up to the trivial (N e + l)-fold degeneracy in 
finite-size clusters with open boundary conditions. 

Before investigating the above properties, we consider the t-ladder at t — 0. The ground 
states of this model can be written in the following form: 

L 

l*) = ®|a)i, (2-1) 

i=l 

where |a)i's correspond to either of states (i)-(iii) defined in table [I] for n < 0.5, and (ii)-(vii) 
for n > 0.5. The degeneracy of the ground states, the energy E and the chemical potential 
/i(= dE/dN e ) are summarized in table [TT]. The charge gap A c at quarter-filling {n = 0.5) is 
2t±. 

Next, we consider the t-ladder in the limit t/t± — > 0. Let us consider the cases n < 0.5 
and n > 0.5, separately. 
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A. n < 0.5 

Up to order t, the effective Hamiltonian is 

Kf = -t^(6i& !+lff +h.c), (2.2) 

icr 

where b\ a denotes a creation operator of an electron in the bonding band at rung i with spin 
a {a =t? I) with the constraint that no rung is doubly occupied, or b\ a = b\ a (l — n\ — nf). 
Here, the creation operator b\ a is defined as b\ a = (c^ + c^)/a/2. This Hamiltonian is 
equivalent to the U = oo Hubbard chain. Thus, the ground states are degenerate with 
respect to the spin degrees of freedom. The charge part of the ground-state wavefunction is 
simply that of the spinless fermion model on a chain. 

Next, we consider the effective Hamiltonian of order t 2 /tj_. Here, we define the local 
Hamiltonian 7"^° C +1(T : 

n^+ia = -t(c£cl +la + clc 2 l+la + h.c). (2.3) 
Letting this local Hamiltonian operate on \a)i (g) \(3)i + \, we obtain the following relations: 

H% +la \B - a) t <g> \Ba) i+1 = ^j=\S)i ® |00) m , 



W&.Jfl<r)< ® \B - a) l+1 = ^|00>, ® \S) i+1 , 



V2 l 
at 

72' 



KiM-Mi ® |00) i+1 = ^(\B- <r)i ® \B<r) i+1 - \A - a), ® |Ar) i+1 ), 
^!° c + i CT |00) i ® = --^(15(7), ® \B - a) i+1 - \Aa)i ®\A- a) i+1 ), 

where \S)i and \Aa)i are defined as 4|(c^c^| — cj|c^)|00)i and -^(c^ — c^)|00)j, respectively. 
Thus, the second-order perturbation energy E 2 is obtained as in table |T|. It is also shown 
that the following relations are satisfied: 

n%n? +u \B t)i ® |B = |B t>< ® |S |) m , 
t)< ® \B = |B |>i ® |B T) m , 

where = bLhi and = blJbia- Thus, the effective Hamiltonian TC eS , up to order t 2 /t±, 
is written as follows: 



1 

? 



^eff(2) = ^^t^ ^^^. ^^^^ _ l\_ la nf_ a b i+la + h.c.) 
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where is the spin operator in the bonding band at rung i and J e g = t 2 /t±. Hence, it is 
shown that the effective Hamiltonian, up to order t 2 /tj_ for n < 0.5, has the same form as 
that of the U — > oo Hubbard chain, up to order t 2 /U for n < 1 fl32"|-|35[ 



ground-state properties of the t-ladder in the limit t±/t 
those of the U 



oo Hubbard chain for n < 1 [35 



As a result, the 
oo for n < 0.5 are the same as 
This leads to a spin-singlet ground state 



by the Lieb-Mattis theorem [36 



B. n > 0.5 



In this subsection, we consider the case n > 0.5. The unperturbed ground states are 
written in the form of eq. (|2.1|) . The matrix elements of the local Hamiltonian eq. ( |2.5|) are 
summarized in table [TV]. 



n 



loc 

M+l 



t^2(cllc l i+la + cf a c 2 +la + h.c). 



(2.5) 



Considering the matrix elements in table it is shown that the state \B — a)i (g |crcr) i+1 can 
reach the state |erer)j eg \B — cr) i+1 after successive multiplication by the local Hamiltonian 
Ti l °i + i as follows: \B — a)i eg |crcr)j +1 — > [a — cr)j eg |£>cr)i + i — >■ |5cr)j ® \a — cr) i+ i — > |crcr)j ® 
[5 — Using this property, we can show that the following processes are possible by 

applying the local Hamiltonians j and ?^J° C +1 successively. 



® |ccr)i <8> | cr — cr 
|-Bc)i-i ® |o"cr)i ® | — oo 
\Bo)i-\ eg) | — cr — o)i ® | — oo 
\Bo)i-\ eg) | — cr — o)i Cg) I cr — cr 
\Bo)i-i eg) |crcr)j eg) I — o — o 
|Scr)j_i eg \o — o)i Cg | — oo 



i+ i «-> l-BoVi (g 
i+ i ^ |-Bcr)i_i Cg) 
i+i «-> |-Bo-)i_i (g 
i+i «-> l-Bo-)^! (g 
i+ i ^ |-Bcr)i_i Cg) 

i+ i «-> |s«7)i_i eg 



a - cr)i (g |crcr) i+1 , 

- oo)i eg |crcr)j + i, 

- atr)i (g | - cr - o) i+l , 
cr - cr)i <g | - cr - cr) i+ i, 

- (J - cr)i eg |crcr) i+ i, 

- crcr)i (g [cr - cr) m . 



Thus, | If), | tl)j I IT) and | ||) can have their positions changed after successive multi- 
plication by the Hamiltonian (Tit), if there exists at least one \B j) or \B {). Furthermore, 
the following processes are also possible: 



\oo)i^ ® \Bo)i <g \B 
\o - er)i_i ® \Bo)i eg \B 



o)i+\ <-> |crcr)j_! eg \B - o)i eg \Bo) i+ i, 
cr)i+i <-+\cr- cr)j_! (g \B - a)j eg \Bo) i+1 . 



Thus, \B |) and \B {) can have their positions changed, if there exists at least one | ft), 
I ID? I IT) or I ID- As a result, any unperturbed ground state in the form of eq. Q2.1| ) 
can be reached after successive multiplication by the Hamiltonian (Tit) for 0.5 < n < 
1 in the subspace of fixed total S z and fixed electron number. This property is usually 
called connectivity. Combined this property and the fact that the matrix elements of the 
Hamiltonian in this representation are all positive (table [TV]) , the Perron- Frobenius theorem 
ensures that the state of the largest eigenvalue is unique and the wavefunction is positive 
(nodeless) in this representation in each subspace. 

Next, we consider the ground-state wavefunction in the subspace of the maximum total 
S z (S z = S max ). The Hamiltonian in this subspace has the same form as that of the spinless 
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fermion model on a chain. Thus, the wavefunction of the largest eigenvalue in this subspace is 
nothing but that of the spinless fermion model (|\&sf))- Applying the spin-lowering operator 
S~ to I^sf), w e obtain the eigenstates of various total S z, s. These states have all positive 
(nodeless) wavefunctions in the present representation. Thus, these states have finite overlap 
with the states of the largest eigenvalue in the corresponding subspaces. By use of the gauge 
transformation cf — > (— l) l c°, a = 1, 2, the sign of the hopping amplitude t can be changed, 
i.e. t — > —t, with spin operators unchanged. As a result, it is shown that the ground state of 
the t-ladder in the limit t±/t — > oo for 0.5 < n < 1 has the maximum total spin (S = S mSuX ) 
and is unique up to the trivial (N e + l)-fold degeneracy. 



C. Remarks 



Here, we give some remarks on the above theorems. 



1. Part of theorem [I] can be extended to higher dimensions. The effective Hamiltonian of 
double-layer i-models up to order t 2 /t± for n < 0.5 has the same form as that of the 
single-layer Hubbard models up to order t 2 /U for n < 1 p2]-p4|. Thus, the ground- 
state properties of double-layer t-models in the limit t±/t — > oo for n < 0.5 are the 
same as those of the single-layer U — > oo Hubbard models for n < 1 [[37 . 



2. The proof of theorem is mathematically similar to that of Kubo's theorem 38 



However, the physical situation is different. In Kubo's theorem, the limit of the strong 
Hund-coupling is taken, i.e. i^Hund = &l ' Jn °°- Furthermore, almost 

degenerate bands are assumed. On the other hand, in theorem |2|, we do not assume 
explicit ferromagnetic couplings. In contrast to Kubo's case, the limit of the large 
band splitting is taken in the present case. The extension of Kubo's theorem to 
n < 0.5 shows that the ground state is also ferromagnetic which is contrasted 



with theorem |I[ The proof of theorem is mathematically similar to that of ref. |40 
for the one- dimensional Kondo-lattice model, too. 

3. The restriction on the boundary condition can be relaxed such that the Hubbard chain 
has a unique spin-singlet ground state for theorem [I] and that the spinless fermion 
model on a chain has all positive matrix elements in the site representation for theorem 
H]. For example, theorem |2| can be extended to the case of periodic boundary conditions 
with an odd number of electrons and an even number of rungs. Nagaoka's theorem is 



recovered for the one-hole case 14-16 



III. FERROMAGNETISM 

In this section, we present the numerical results on the t-ladder for finite t±/t and the 
t-J ladder for small J and Jj_ obtained by the DMRG method (finite-size algorithm) ||25|| . 
The DMRG calculation has been performed with open boundary conditions. 

As a test of the DMRG calculation, we compare the ground-state energy obtained by 
the DMRG method with that of the exact diagonalization method. In Fig.|I|, the agreement 
of the data obtained by these two methods is good. The maximum error is about 0.01%. 
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Next, we consider the truncation error, i.e. the error due to small m, where m is the number 



of states kept in the superblock f25j . The difference between m = 50 and m = 100 is very 
small as shown in Fig.|2|, suggesting that m = 50 is sufficient. (See also Fig||.) Thus, we 
mainly report the results of m = 50 hereafter. 

Let us consider ferromagnetism of the t-ladder. In Fig. |3|, we show the energy difference 
AE-p between the ground-state energy in the subspace of S z = and that of S z = S max as a 
function of filling. The data in various-size clusters are scaled to a single line, indicating that 
the finite-size error is small. (See also Fig.[7|.) From this plot, we find the region FF where 
the fully-polarized ferromagnetic state is one of the ground states. The phase boundary n ci 
between FF and non-FF is estimated as shown in Fig. |TT. At t±/t — 1, the result in ref. f2H 



(n Cl ~ 0.78) is recovered. Qualitatively, the region FF becomes larger as t±/t increases. The 
phase boundary n Cl gets closer to 0.5 as t±/t increases. This is consistent with the rigorous 
results in Sec.|l[ 

Next, we consider the phase boundary n C2 between SS and non-SS, where SS is defined 
as the region in which the ground state is a spin-singlet. Figure f| shows the ground-state 
energy as a function of the total S z at t±/t = 1 for n = 0.5625, 0.625 and 0.6875. The ground 
state is a spin-singlet for n = 0.5625 and not for n = 0.625 and 0.6875. Hence, the phase 
boundary n C2 is estimated as n C2 = 0.59 ± 0.03, which is consistent with ref. |M](n C2 ~ 0.6). 
In this way, the phase boundary n C2 is estimated for various t±/t as shown in Fig.[TT]. The 
region PF shrinks as t±/t increases, where PF is defined as the region which is neither FF 
nor SS. This is also consistent with the rigorous results in Sec.[n[ 

Here, we present the numerical results on the t-J ladder for small J and J±. For sim- 
plicity, we choose t± = t and J± = J, and set t = 1 as the energy unit. Figure || shows the 
n-dependence of AE F at J = 0.00,0.05,0.07,0.10 and 0.15. As shown in this figure, AE F 
becomes large near half- filling (n = 1) due to antiferromagnetic correlation, and AE-p seems 
to be the smallest near n ~ 0.8 for n > 0.5. Figure |5| also indicates that J = 0.05 is enough 
to destroy the region FF. 

Next, we consider the stability of PF against J. The ground-state energy as a function 
of the total S z near half-filling for J = 0.05,0.07,0.10 and 0.15 is shown in Fig. || This 
figure suggests that FF is surrounded by PF in the phase diagram of the t-J ladder for finite 
5. 



IV. METAL-INSULATOR TRANSITION AT QUARTER-FILLING 

Before discussing the metal-insulator transition, we consider the charge gap A c at 
quarter-filling (n = 0.5). Figure [7| shows the n-dependence of the chemical potential u 
(= 8E/dN e ). The chemical potential u in a finite-size cluster is defined as 

s EM -EM 
[ni -n 2 )N s 

where E(rii) denotes the ground-state energy at filling rii, i = 1,2, and n = (rii + n 2 )/2. 
We took (n 2 — ni)N s = 2. The charge gap A c is defined as A c = u(n c + 0) — u(n c — 0), 
where n c is the critical electron-density (n c = 0.5 in the present case). It is expected that 
the charge gap opens at quarter-filling, if t±/t is large enough (Sec.||). Actually, for large 
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values of t±/t, the charge gap seems to open as shown in Fig.|7|(a). For smaller values of 
t±/t, we cannot determine whether the charge gap opens from Fig.[7[ Thus, we extrapolate 
the charge gap in a finite-size cluster |eq.( |4.2| )l as a + b/L, using the data for L = 12 — 24, 
and estimate the charge gap as shown in Fig. || 

A C (N, = JV./2) = E(N ' = N ' /2 + 2) + E(N " = 2 NJ2 ~ 2) " 2E(N ' = N ' l2) . (4.2) 

There are some possibilities for the critical value t± c /t above which the charge gap 
opens. One of the possibilities is that the critical value t± c /t is zero and that the gap is 
exponentially small in the limit t±/t — > as in the case of the Hubbard chain in the limit 
U — > pBfl . In order to determine the critical value t± c /t, we adopt the extended Aharonov- 
Bohm ( AB) method proposed by Kusakabe and Aoki . In the framework of this method, 
we investigate the extended spectral flow by introducing a Peierls phase as 

C^+i. - expO^cga^. (4.3) 

It is expected that the period of the spectral flow reduces to L^q/M if M-particle bound 
states are formed in the ground state. For example, M = 1 for a metallic state, and M = 2 
for a BCS state. We apply this method to the t- ladder with a very small value of t±/t 
(t±/t = 0.001). As shown in Fig.|9], the spectral flow at quarter-filling has the minimum 
extended AB period, i.e. $ = $ , suggesting that the ground state is an L-particle bound 



state or, in this case, an insulator [42]. This behavior is contrasted with the case off quarter- 



filling. For example, at n = 1/3, the extended AB period is larger than $ as shown in FigJlO. 
This implies that the ground state at quarter-filling is an insulator already for t±/t = 0.001. 
It is plausible to consider that the critical value t± c /t may be zero. A possible scenario may 
be the following: The perturbation of the small tj_-term produces the relevant Umklapp 
process which leads to an insulator, at the same time as the degeneracy with respect to 
the spin degrees of freedom is removed. The numerical results presented above are quite 
different from the case in the weak-coupling limit (U — > +0). In the weak-coupling limit, 
it is shown by bosonization that t± c /t is one J7|. Thus, it is expected that t± c /t decreases 
from 1 to as the interaction U increases from to oo. 

Now, let us consider the metal-insulator (MI) transition at quarter-filling. As discussed 
in Sec. O, in the limit t±/t — > oo, the MI transition for n — ► 0.5 — is effectively described by 
the equivalent model to the U — > oo Hubbard chain, and the MI transition for n — > 0.5 + 
is described by the spinless fermion model on a chain. It is interesting to compare these 
features with those in the weak-coupling regime. In the weak-coupling regime (U — > +0), the 
charge gap is also expected at quarter- filling for t±/t > 1 because of the relevant Umklapp 
process . This MI transition is understood as the Mott transition which is described by the 
U — > +0 Hubbard model on a chain written in terms of the bonding-band operators. In both 
weak-coupling [U <^ t±(> £)] and strong-coupling [U 3> t±(^> t)] regimes, as n — >• 0.5 — 0, 
the MI transition is described by single chain effective Hubbard Hamiltonians. However, 
the value of the charge gap will have different energy scales. In the weak-coupling regime, 
the value of the charge gap would be determined mainly by U. On the other hand, in the 
strong-coupling regime, it would be determined mainly by t±. This feature is similar to the 
two types of the MI transition for transition-metal compounds |43[], i.e. the Mott-Hubbard 
type and the charge-transfer type. 



S 



V. SUMMARY 



In summary, two aspects of the ground state of the U — oo Hubbard ladder are inves- 
tigated. One is ferromagnetism, the other is the metal-insulator (MI) transition. In the 
limit t±/t — > oo, it is rigorously shown that the ground state is a spin-singlet for n < 0.5 
and that the total spin is maximum for 0.5 < n < 1. For finite t±/t, we have estimated 
the phase boundaries, with respect to spontaneous magnetization, by the density-matrix 
renormalization group method. It is numerically shown that the region FF becomes larger 
and spreads down to quarter-filling as t±/t increases, which is consistent with the rigorous 
results presented in Sec.|TJ. The rigorous results (t±/t — > oo) and the numerical results for 
finite t±/t support one another and confirm that the ground state can be ferromagnetic 
for the U = oo Hubbard ladder with finite hole-density. The numerical results for the t-J 
ladder suggest that FF is surrounded by PF for finite 5 in the small J regime. We have also 
estimated the value of the charge gap at quarter-filling (n = 0.5). Applying the extended 
Aharonov-Bohm method, we have obtained numerical evidence that the critical value t± c / 1, 
above which the charge gap opens, is less than 0.001. This is quite different from that of 
the weak-coupling limit (U — > +0) (t± c /t = 1) 0. 
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TABLES 



TABLE I. Basis set. 





SvvnHnl 

y ±l±uvj 1. 


T~)pfi n it, inn 




(i) 


\00}i 


vacuum 





(ii) 




^(cjt + c^loo), 


-*± 


(iii) 


\Bl)i 


72(41 + ^)100), 




(iv) 


1 n>i 


cJfc^OO), 





(v) 


IU)i 


cJcJlOO), 





(vi) 


\nh 


£?5j|00)i 





(vii) 


lit), 


cJ4 f |00), 






TABLE II. Ground states at t = 0. 





n < 0.5 


0.5 < n < 1 


Degeneracy 
Energy £7 
Chemical Potential /i 


L C Ne x 2^ 
-*± x N e 
-t± 


L CAr e _ L x 2^ 
-t_L x (2L - iV e ) 
t± 




TABLE III. Second-order perturbation enerj 




£?2 between \a) and /3 


) |a> 


1/3) 




-t 2 /(2t ± ) 
t 2 /(2t ± ) 
t 2 /(it ± ) 

-ty(4t ± ) 


|5<7)i 
|B<7)j (8) \B- a) i+ i 
\Bo)i®\B-o) l+1 
|00)i_i |S<7)i ® |B - 
|00)i_i |Bff)i |B - <7) i+ i 


® |5<7) i+ i 
|£?cj)i <g> \B - a) i+1 
\B-a)i <g> \Bo) i+1 
\Bo)i ® \B - a)i <g> |00) i+ i 
\B - a)i ® |5<7)i <g> |00) i+ i 
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TABLE IV. Matrix elements. 



/ ll_/loC 1 Q\ 


\a) 


1 /3\ 
\P) 


t 


\Ba)i (g> |<Jcr)j + i 


|(7cr)j (g) |-B(7) i+ i 




\Bo)j ®\a — cr),-4-i 


1(7 — 0"),- (8) |i?(j),--L-| 




\Ba)i <g> |cr - cr) i+ i 


\aa)i ®\B- a) i+1 


t/2 


SCT^ (g) - (7C7)i + l 


- aa)i <g) \Ba) i+ i 




-B(j)i <g) - (7C7)j + l 


\aa)i (g> \B - a) i+ i 




® | - o - a) i+ i 


- aa)i (g> \B - o) i+ i 




\Bo)i <g> | - a - a) i+ i 


\a - a)i (g> \B - a) i+ i 





otherwise 
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FIGURES 

FIG. 1. Energy per site as a function of the total S z , measured from that of S z = S max , for the 
t-J ladder at J = J± = 0.05 and t = t± = 1 in a 10 x 2-site cluster with 18 electrons. The solid 
diamonds and open squares correspond to the data obtained by the DMRG method and the exact 
diagonalization method, respectively. For the DMRG method, we took m = 50 (m : the number 



of states kept in the superblock |p5|) and repeated 2-3 sweeps for convergence. 



FIG. 2. Energy difference AE-p per site between the ground-state energy in the subspace of 
S z = and that of S z = S max as a function of filling for the t-ladder at t±/t = 1. The solid 
diamonds and open squares correspond to the data of m = 50 and m = 100, respectively, t = 1. 

FIG. 3. Energy difference AE-p per site as a function of filling for the t-ladder at t±/t =2.5 
[(a)], 1.0 [(b)] and 0.5 [(c)]. t = 1. 

FIG. 4. Energy per site as a function of the total S z , measured from that of S z = S max , for 
the t-ladder at t = t± = 1 in a 16 x 2-site cluster with 22, 20 and 18 electrons starting from above. 
The solid and open symbols denote the data for m = 50 and m = 100, respectively. 

FIG. 5. Energy difference AEp per site as a function of filling for the t-J ladder in 12 x 2-site and 
16 x 2-site clusters at J/t = J±/t = 0.00, 0.05, 0.07, 0.10 and 0.15 starting from above, t = t± = 1. 

FIG. 6. Energy per site as a function of the total S z , measured from that of S z = S max , for 
the t-J ladder at t = t± = 1 in a 16 x 2-site cluster with 26 [(a)], 28 [(b)] and 30 [(c)] electrons 
at J/t = J±/t =0.05, 0.07, 0.10 and 0.15 starting from above. The dotted lines correspond to the 
energy of S z = 0. The data are normalized by J. 

FIG. 7. Chemical potential /i as a function of filling at t±/t =2.5 [(a)], 1.0 [(b)] and 0.5 [(c)]. 
For comparison, the chemical potential [i for the non-interacting case (dotted line) and that of the 
spinless fermion model (solid line) are shown {N s = 160 x 2). The points A and B correspond to the 
anomalies due to the band bottom of the anti-bonding band and the band top of the bonding band, 
respectively. The points F and S correspond to the phase boundaries n Cl and n C2 , respectively, 
which are estimated by the DMRG method (Sec.pII|). t = 1. 



FIG. 8. (a) Size dependence of the charge gap A c for the t-ladder at t±/t =2.5, 1.5, 1.2, 1.0, 
0.8, 0.7, 0.5 and 0.2 starting from above, (b) The charge gap A c as a function of t±/t. The bold 
line is a guide to the eye. t = 1. 

FIG. 9. Spectral flow of the i-ladder for t±/t = 0.001 in a 12-site cluster with 6 electrons for 
< <I> / ^ 1 [( a )]> an d the blow-up region for <I> / — 0.5 [(b)]. The solid diamonds correspond to 
the spectral flow of the ground state. The data are obtained by the exact diagonalization method. 
t = 1. 
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FIG. 10. Spectral flow of the t-ladder for t±/t = 0.001 in a 12-site cluster with 4 electrons for 
< ^/^o < 1 [( a )]> an d the blow-up region for <I> / — [(b)]. The solid diamonds correspond to 
the spectral flow of the ground state. The data are obtained by the exact diagonalization method. 
t = 1. 

FIG. 11. Phase diagram of the U = oo Hubbard ladder with respect to spontaneous magne- 
tization. The dashed lines denote the region where the ground state is rigorously shown to have 
the maximum total spin S (S = S ms , x ) (Theorem [2] and Nagaoka's theorem). The dotted line in 
the limit t±/t — > oo corresponds to the region where the ground state is shown to be a spin-singlet 
(Theorem [l]). The solid and open diamonds correspond to the phase boundaries n Cl and n C2 , es- 
timated by the DMRG method. The bold lines are guides to the eye. At t±/t = 0, the ground 
states are degenerate with respect to the spin degrees of freedom, because the model reduces to 
the decoupled U = oo Hubbard chains. At quarter-filling, the charge gap is expected (Sec,[TV|). 
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